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Abstract. Practical expressions are derived for evalua- 
tion of electrical and thermal conductivities and ther- 
mopower of degenerate electrons in the outer envelopes of 
neutron stars with magnetic fields. All tensor components 
of the kinetic coefficients are calculated (those related to 
conduction along and across magnetic field and to the Hall 
currents). The kinetic coefficients are presented as energy 
averages of expressions containing energy dependent effec- 
tive relaxation times of two types, associated either with 
longitudinal or with transverse currents. The calculation 
is based on the effective scattering potential proposed in 
the previous paper, which describes the electron-ion and 
electron-phonon scattering, taking into account correla- 
tion effects in strongly coupled Coulomb liquid and multi- 
phonon scattering in Coulomb crystal, respectively. An- 
alytic fitting formulae are devised for the effective relax- 
ation times at arbitrary field strength. Basing on these 
results, we calculate the transport coefficients at various 
temperatures, densities, and magnetic fields pertinent to 
the neutron star envelopes. 
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1. Introduction 

Transport properties of neutron star envelopes determine 
many aspects of neutron-star evolution. For instance, the 
thermal conductivity in the outer envelope affects cool- 
ing of a neutron star and its radiation spectra (e.g., Gud- 
mundsson et al. |1983j; Page |1997|; Potekhin et al. |1997|). 



The electrical conductivity is the basic quantity for the 
studi es of magnetic-field ev olutio n (e.g., Muslimov & Page 
1996; Urpin & Konenkov 1997; Konar & Bhattacharya 



1997), w hich in turn affects thermal evolution (Miralles et 



al. 1998 ). The thermopower determines a v ariety o f ther- 
momagnctic phenomena (Urpin & Yakovlev 1980b; Urpin 



et al. [1986| Shibazaki et al. |l989| ; Yabe et al. |l99l| ). 

In the outer envelopes of neutron stars, the trans- 
port coefficients are mainly determined by the processes 
of electron scattering off strongly correlated ions. Gen- 



eral formalism for calculating kinetic properties of strongly 
coupled Coulomb plasmas was developed by Hubbard & 



Lampe (1969) and Flowers & Itoh (1976) (see references to 
earlier results therein). Yakovlev & Urpin ( 1980] ) derived 
approximate analytic expressions, which were confirmed 
later in more elaborate calculations by Raikh & Yakovlev 
(1 19821 ), Itoh et al. (|1983| ) and Nandkumar & Pethick 
(1984). Itoh et al. Ql984[ |1993|) improved the results of 



Yakovlev & Urpin ( p8q ) and Raikh & Yakovlev ( p82| ) 
in the solid crust by taking into account finite sizes of 
atomic nuclei (which may be important in the inner crust) 
and the Debye- Waller factor (which describes reduction of 
electron-phonon scattering rate due to background lattice 
vibrations). The Debye- Waller factor proved to be impor- 
tant at temperature T close to the melting temperature 
of a Coulomb crystal, T m , or at sufficiently high densities 
where zero-point lattice vibrations are strong. Detailed 
numerical and analytic calculations by Baiko & Yakovlev 



( |1984j [1993D 



1995 , I99C) were in reasonable agreement with Itoh et al. 



Magnetic fields in the neutron-star envelopes com- 
plicate electron transport making it, particularly, anis- 
otropic. The field strengths of radio pulsars range from 
B ~ 10 8 G to B > 10 13 G, with typical value B ~ 10 12 G 
(Taylor et al. 1993). Some X-ray pulsars and soft gamma 
repeaters are probably magnetars - neutron stars with 
B ~ 10 14 - 10 15 G, as suggested by Thompson & Dun- 



can (1995) and supported by recent observations (Va- 
sisht & Gotthelf 119971; Gotthelf et al. 119991; Kouveliotou 



19981 , [199E 



Shitov 1999). Magnetic field strength 



et al. 

expressed in the relativistic units, 
b = heB/(m 2 c c 3 ) w S/(4.414 x 10 13 G), 



(1) 



is greater than unity for magnetars, unlike for ordinary 
pulsars. 

The magnetic field affects thermodynamic and kinetic 
properties of dense degenerate plasmas in different ways, 
depending on density, temperature, and field strength 



(e.g., Yakovlev & Kaminker 1994). In general, electron 



motion transverse to the field is quantized into Landau 
orbitals. For sufficiently high temperature, however, the 
field can be treated as non- quantizing (classical). The 
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non-quantizing magnetic field does not affect thermody- 
namic properties of matter, but hampers transverse trans- 
port and causes the Hall currents. In a weakly quantiz- 
ing field, where electrons populate several Landau levels, 
the thermodynamic functions and kinetic coefficients oscil- 
late with increasing density around their classical values. 
Strongly quantizing magnetic field confines most electrons 
to the ground Landau level. In this case, thermodynamic 
and kinetic properties of matter are very different from 
those in the classical regime. 

The problem was extensively studied since the end of 
1960s (Canuto & Ventura 1977 and references therein). 
Expressions for the kinetic coefficients, based on the so- 
lution of the Boltzmann equation in the relaxation-time 
approximation and valid for non-relativistic and relativis- 
tic electrons in the neutron star envelopes with quantizing 
magnetic fields were obtained by Kaminker & Yakovlev 
( |l98lD for transverse transport, by Yakovlev (1984) for 
longitudinal transport and by Hernquist (1984) for all 
cases. Potekhin (1996, hereafter Paper I) has shown that 
the usual description of electron transport along magnetic 
field with the Boltzmann equation yields almost the same 
transport coefficients as a more general kinetic equation 
for the electron spin density matrix. According to Paper I, 
the Debye- Waller factor enhances the longitudinal electri- 
cal and thermal conductivities near T m much stronger if 
the magnetic field is q uantiz ing. Basing on these results, 
Potekhin & Yakovlev ( 1996 , hereafter Paper II) derived 
practical formulae for the longitudinal transport coeffi- 
cients. 

The studies cited above used the customary one- 
phonon approximation for the electron scattering in 
Coulomb solid. Recently Baiko et al. (1998) have recon- 



sidered electron transport in non-magnetized plasmas by 
including multi-phonon processes, which have proved to 
give a contribution of similar magnitude (but opposite 
sign) as the Debye- Waller factor near T m . Concerning the 



Coulomb liquid, Baiko et al. (1998) have suggested that 



incipient ordering of ions in the strong-coupling regime 
affects electron scattering; they have proposed an approx- 
imate treatment of this effect by modification of the static 
structure factor of ions. Both modifications (in the solid 
and liquid phases) change the kinetic coefficients near the 
melting point and drastically reduce their discontinuities 
at T — T m . The new ion structure factors have been em- 
ployed by Potekhin et al. (1999, hereafter Paper III) in 



calculations of electrical and thermal conductivities in the 
outer envelopes of neutron stars without magnetic fields. 
Numerical results in Paper III have been fitted by analytic 
expressions, derived in the relaxation-time approximation 
with the use of a specially adjusted effective scattering 
potential. 

In this paper, the effective potential obtained in Pa- 
per III is applied to calculation of electron transport co- 
efficients at arbitrary magnetic field strength, at tem- 
peratures T ~ (10 5 - 10 9 ) K and densities p ~ (10 3 - 



10 11 ) gcm~ 3 typical for the outer envelopes of neutron 
stars. Energy-dependent effective relaxation times, sub- 
ject to thermal averaging, are obtained for electron trans- 
port parallel and perpendicular to the quantizing mag- 
netic fields. In the case of non-quantizing fields, the usual 
semiclassical formulae (e.g., Urpin & Yakovlev 1980b) are 



utilized taking into account the results of Paper III. 

The paper is composed as follows. In Sect. ||, we de- 
scribe typical plasma parameters of interest. In Sect. [| we 
express the electron transport coefficients through an ef- 
fective scattering potential in the relaxation time approxi- 
mation. In Sect. [|, we present analytic fits to the effective 
energy-dependent relaxation times related to longitudinal 
and transverse electron transport in quantizing magnetic 
fields. Section || illustrates the main features of the trans- 
port coefficients given by the present theory. 



2. Magnetized degenerate matter 

Consider a plasma composed of electrons (with charge — e) 
and a single ion species with charge Ze and mass m\ w 
Am u (where m u = 1.6605 x 10~ 24 g is the atomic mass 
unit, and A is the atomic mass number). The complete 
pressure ionization occurs at high temperatures or high 
densities (e.g., at p J; 22 Z 2 A gcnr 3 in the non-magnetic 
case - see Potekhin et al. 1997). Whenever this assumption 
is violated, we will employ the mean-ion approximation, 
in which all ions in all ionization stages are replaced by 
a single species with some effective values of Z and A. 
Electrons are assumed to be degenerate and nearly free. 
The degeneracy implies T < Tp = (ep — m c c 2 )//cb, where 
ep is the Fermi energy (including the rest energy, m c c 2 ) 
and fee is the Boltzmann constant. Degenerate electrons 
can be considered as nearly free, if their kinetic energy 
exceeds greatly a typical energy of electron-ion Coulomb 
attraction; in a non-magnetized plasma, this happ ens at 
p > 10AZ gem' 3 (e.g., Pethick & Ravenhall [l995| ). 

A degenerate electron gas can be characterized by the 
Fermi momentum pp or wave number /cp = pp/h. With- 
out any magnetic field, kp = /cpo = (Sn 2 ^) 1 ^ , where 
n c pZ I (m u A) is the electron number density. It is also 
convenient to introduce the "relativity parameter" (e.g., 
Shapiro & Teukolsky |1983| ) 



x T = hk-po/m c c ~ 1.009 (p^Z/A) 



1/3 



(2) 



where pe = p/(W 6 gcm~ 3 ). 

Let the magnetic field B be directed along the z-axis. 
Then, using the Landau gauge of the vector potential, the 
quantum states of a free electron can be labelled by the 
y-coordinate of the electron guiding centre, the longitu- 
dinal electron momentum p z , the Landau number n, and 
a spin variable s. The ground Landau level (n = 0) is 
non-degenerate with respect to the spin variable (s = — 1, 
statistical weight go = 1) while the other levels (n > 0) 
are doubly degenerate (s = ±1, g n = 2). 
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It is convenient to write (Paper II) 



where e is the electron energy, 
/o(e) = {exp[(e- / i)/fc B T] + l}- 1 



(3) 



(4) 



is the Fermi-Dirac distribution, p is the chemical potential 
(e and \i include m e c 2 ), and 



(5) 



Here and hereafter, w c = eB/m D c is the electron cyclotron 
frequency, p n (e) is the value of \p z \ for an electron on 
the Landau level n, and n max is the maximum Landau 
number for a given energy e. Since e = (m 2 c 4 + c 2 p 2 + 
2m e c 2 %w c n) 1 / 2 , we have 



Pn(e) = [(e/c) 2 - (m c c) 2 - 2m e hoj c n] 1 ^ 2 



(6) 



and obtain n max as integer part of convenient (Paper I) 
dimcnsionless energy variable 



v = Po( e )/( 2m e^ w c)- 



(7) 



At T <C Tp, the derivative (—dfo/de) can be replaced 
by the delta function 5(e — cf); in this case n c = A/b(ef)- 
If "max ^ 1 , the sum in Eq. (||) can be approximated by 
an integral, which gives the classical (field-free) result, 



A/o(e)=^(e)/(3rf). 



(8) 



Electrostatic screening produced by electrons is char- 
acterized by the Thomas- Fermi wave number /ctf^ 



47re' 



dn c 
dp 



47re z 



<9e 



<9e 



de. 



(9) 



The field can be considered as non-quantizing if T 
Tb , where (Yakovlev & Kaminker 1994) 



T B = huj g /k B » 1.34 x 10 8 (B 12 / 7r ) K, 



(10) 



B 12 = B/(10 12 G), 7 r = x/l + a;?, and w g = eSc/e is the 
electron gyrofrequency (for e = £f). In the non-quantizing 
fields, T F = T r ( 7r - 1) and (fc TF /2fc F o) 2 « (a f /tt) 7r /x r , 
where T r = m c c 2 /fc B « 5.93 x 10 9 K and a f = e 2 /Tic w 
1/137 is the fine-structure constant. 

The opposite case of strongly quantizing field occurs 
at T <C Tg and p < pb, where 

Pb = Am u n B /Z « 7.045 x 10 3 (A/Z)S 3 ^ 2 gem" 3 , (11) 

tib = (ttv^) -1 s m 3 , and a m = {Uc/eB) 1 / 2 is the so called 
magnetic quantum length. In this case, 

fc F - n c = (4/3) 1 / 3 ( j0 /p B ) 2 / 3 k F0 . (12) 

Therefore Tp is strongly reduced for p <C Pb ■ 



The state of the one-component plasma (OCP) of ions 
depends on the Coulomb parameter, 



(Zef 



22.75 Z 2 

Tr 



1/3 



(13) 



where a\ = [3/(47rn i )] 1 / 3 is the ion-sphere radius, n\ = 
n e /Z is the number density of ions, and T 6 = T/(10 6 K). 
In a weakly coupled OCP, T <C 1, ions form the Boltzmann 
gas whose screening properties are characterized by the 
inverse Debye screening length, 



(14) 



For L > 1, the ions constitute a strongly coupled liquid. 
The liquid freezes into a Coulomb crystal at some Y = Y m . 
For classical ions (whose zero-point quantum vibrations 
are negligible), r m w 175, whereas strong zero-point vi- 
brations suppress the freezing and increase r m (Nagara et 
al. 1987 ). The freezing is completely suppressed in the so 
called quantum liquids, which exist at x T 0.18AZ 7 / 3 , 
as can be estimated from numerical simulations (Jones & 
Ceperley 1996). In general, the quantization of ionic mo- 
tion is significant at T <C T p , where 



T p = Tuv p /k B w 7.832 x 10 b (Z/A) ^Pe K 



(15) 



is the ion plasma temperature, w p = ( / ±irZ 2 e 2 n\/m/ l ) 1 / 2 
being the ion plasma frequency. We do not consider the 
quantum ion solids and liquids hereafter. 

We neglect also effects of magnetic field on the OCP 
of ions. This is justified if the ion cyclotron energy huj C i = 
hoj c Zm c /mi is small compared with either k B T [i.e., Tq ^> 
0.0737 (Z/A) Bi 2 ] or typical phonon energies in the OCP 
(~ hto p ; i.e., y/pt > 0.0094 B 12 ). 

The characteristic p — T domains are shown in Fig. [l] 
for iron plasma at B = and 10 12 G. We have taken into 
account partial ionization in the mean-ion approximation. 
Electrons are degenerate below Tp; ions are classical above 
T p ; thus, strictly speaking, our consideration is valid in the 
stripe between the solid lines. In practice it may be reason- 
ably accurate outside this stripe because it incorporates 
thermal averaging and because the quantum effects are 
actually small as long as T <; 0.1 T p . 

The short-dashed contours indicate the region of par- 
tial ionization: the upper contour corresponds to the ef- 
fective charge Z — 20 and the lower one to Z = 15. 
We have evaluated Z in the same manner as Potekhin 



et al. (1997), requiring the equation of state (EOS) of a 
plasma composed of free electrons and ions with the effec- 
tive charge Z to reproduce a "standard" EOS that takes 
into account partial ionization. In the non-magnetic case, 



the OPAL EOS (Rogers et al. 1996) has been adopted as 



such a standard and, whenever necessary, interpolated as 



explained in Potekhin et al. (1997). In the magnetic case, 



we have used the finit e-tem perature Thomas-Fermi EOS 
by Thorolfsson et al. fll99S|) . 
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Fig. 1. Characteristic plasma domains on the p — T plane 
for iron. Upper panel: non-magnetic plasma; lower panel: 
B = 10 12 G. Solid lines show Tp and T p vs. p; upper and 
lower dot-dashed lines correspond to T = 1 and T = 175, 
respectively; short-dashed lines indicate the domains of 
partial ionization. Dotted lines on the upper panel show 
temperature profiles in the envelope of a "canonical" cool- 
ing neutron star (see text) for two values of the effective 
surface temperature, 2 x 10 5 and 2 x 10 6 K. Long-dashed 
lines on the lower panel show Tb and ps and separate the 
regions of strong and weak magnetic quantization. 



The dotted curves on the upper panel (marked "NS 



"NS") 
1997l )Fl 



reproduce the temperature profiles (Potekhin et al. 199 
in the envelope of a "canonical" neutron star of the mass 
1.4 Mq and radius 10 km, with an effective surface temper- 
ature 2 x 10 5 K (the lower curve) and 2 x 10 6 K (the upper 
curve). In reality, effective temperatures of the middle- 
aged isolated neutron stars are believed to lie between 



these two extremes (e.g., Page 1998) 



We have recalculated the profiles using the updated non- 
magnetic conductivities (Paper III), but an effect of the update 
turned out to be negligible. 



The dot-dashed lines on both panels correspond to T = 
1 (gas/liquid smooth transition, upper lines) and T = 175 
(liquid/solid phase transition, lower lines). 

Finally, the long-dashed lines on the lower panel in- 
dicate three p — T regions, where the magnetic field is 
strongly quantizing (to the left of ps and considerably 
below Tb), classical (much above Tb), or weakly quantiz- 
ing. 

3. Transport coefficients 

3.1. General relations 

Consider electron electric and thermal currents induced in 
a magnetized plasma under the effect of an electric field E 
and weak gradients of chemical potential V/i and temper- 
ature VT. These currents can be decomposed into con- 
duction and magnetization components (e.g., Hernquist 
1984 ). The latter ones relate to surface effects and must 
be subtracted. Let j e and jt be the conduction compo- 
nents of electric and thermal current densities. They can 
be written as (e.g., Landau & Lifshitz 1960) 



3, 



a 



VT, j T 



a ■ 



K 



VT, 



(16) 



E + V/i/e is the electrochemical field. The 
, a, and k denote second-rank tensors (a is 



where E* = 
symbols a, a 

the electrical conductivity tensor) which reduce to scalars 
at B = only. From the Onsager symmetry relation one 
obtains: &ij(B) = Ta oi {-B) = Ta tj (B). 

1 61) can be rewritten as 



E 



Equations 

= R ■ je - Q ■ VT, 



3t = -TQ -j e -K- VT, (17) 



where R — a^ 1 , Q = —R ■ a, and k = k + Ta ■ Q are the 
specific resistance, thermopower, and thermal conductiv- 
ity tensors, respectively (here Q is defined as in Paper II; 
an opposite sign has been adopted by Hernquist 1984 ). 
The components of a, a, and k can be expressed as 







atij 









e(p-e)/T 



AW , J df 



de. (18) 



The functions Ty-(e) play role of relaxation times deter- 
mined by electron scattering in the magnetic field. Owing 
to the symmetry properties of the tensors a, a, and k, 
there are only three different non-zero components: t zz re- 
lated to longitudinal currents, t xx = r yy related to trans- 
verse currents, and r xy — —r yx related to the Hall cur- 
rents. 

3.2. N on- quantizing magnetic field 

If the quantizing nature of the magnetic field is neglected, 
then (e.g., Urpin & Yakovlev [T980b{ ) 



To, 



1 + Kro) 5 



l + (^To) 2 



(19) 



A.Y. Potekhin: Electron conduction in magnetized neutron star envelopes 



791 



where To is the non-magnetic relaxation time, equal to an 
inverse effective collision frequency in this case. 

In the outer envelopes of neutron stars, relaxation is 
mainly determined by electron-ion scattering. We restrict 
ourselves to consideration of this mechanism; possible in- 
clusion of other processes is discussed briefly in Sect. |[ In 
strongly coupled Coulomb plasmas, the scattering is sig- 
nificantly affected by ion correlations. In the liquid phase, 
an appropriate structure factor of ions should be em- 
ployed. In the solid phase, an adequate description is pro- 
vided by the formalism of electron scattering off phonons 
with allowance for multi-phonon processes (see Baiko et 
al. 1998 for discussion and references). In both the liquid 



and solid regimes, it is convenient to write the squared 
Fourier transform of the scattering potential as 



\U„ 



2\2| 



(20) 



Here, <j> q is the so called screening function (which would 
be equal to q~ 2 for the Coulomb potential, were the 
screening neglected). Then 



ro(e) 



PqVq 



47rniZ 2 e 4 A (e) ! 



(21) 



where vo = poc 2 /e is an electron velocity, and Aq(e) is 
the non-magnetic Coulomb logarithm. In strongly cou- 
pled, degenerate Coulomb plasmas, one has 



Ac 



2po/h 



&qq Z \4> q \ 2 S{q)\F{q)\ 2 



hcq\' 



, (22) 



where F(q) is the form factor of ions, and S(q) is an ef- 
fective structure factor that describes the effects of ion 



correlations (Baiko et al. |1998| ) . 

In Paper III two forms of the effective Coulomb loga- 
rithm (An- and A K ) have been obtained from calculations of 
a and k at B — beyond the relaxation-time approxima- 
tion. The latter approximation fails if inelastic processes 
with energy transfe r > k ^T are important. According to 
Baiko & Yakovlev ( 1995 ) and Paper III, this happens at 
low temperatures T <C T p . In this paper we will focus on 
the case T <; T p , in which we can neglect the difference 
between A K and A CT and adopt Ao = A CT for B = 0. 

The transport coefficients due to electron-ion and 
electron-phonon scatterings (in the Coulomb liquid and 
solid, respectively) can be described with a reasonable 
accuracy (Paper III) in a unified way under the condi- 
tions typical for the outer envelopes of neutron stars. For 
this purpose it is sufficient to replace \4> q \ 2 S(q)\F(q)\ 2 in 
Eq. (|22]) by the expression: 



(9 2 + 9 s 2 ) 



w(q) 

-2\2 ^VP'^)j 



where 



£ = {q 2 + k 2 TF )e-^ 



(23) 



(24) 



(1 + 0.06 He"^, 
w(q) =u_ 2 (< Z /g D ) 2 (l + /3/3), 
1 + 0.122/? 2 



G{t v ,p) 



D(t p ), t 



0.19 



(25) 
(26) 

(27) 



In this case e - ™^ 9 ^ plays role of an effective Debye- Waller 
factor at large T and is negligible at T ^ 1; u_2 = 13, 
where Uj is (aj/uj p y averaged over phonon frequencies ui 
in a Coulomb crystal (e.g., Pollock & Hansen 1973; ); q s is 
an effective screening wave number; and G(t p , j3) is a phe- 
nomenological factor that describes reduction of the scat- 
tering rate, caused by quantum effects at t p = T/T p <C 
1. The factors G, w, and q s contain also phenomeno- 
logical corrections to the Born approximation expressed 
through /3 = TrctfZppc/eF. Finally, the function D(t p ) = 
exp[— otoU-i exp(— 9.1t p )/4], where ao = 1.683-^/ x T /(AZ) 
and U-i — 2.8, is associated with quantum corrections to 
the Debye- Waller factor (Baiko & Yakovlev |l995| ). Note 
that one can safely set G = 1 for T ^ T p and Z ^ 30. 

3.3. Transport along quantizing magnetic field 

Let us calculate the longitudinal electron transport coeffi- 
cients in the quantizing magnetic field using the relaxation 
time approximation and the effective scattering potential 
determined by Eq. J23]). According to Papers I and II, the 
longitudinal kinetic coefficients can be written in the form 
( |l8| ) by defining the effective relaxation time t zz = t\\ as 



A/"s(e 



(eB) 



4TT 3 hZ 2 e 4 ni 



He)- 



(28) 



The dimensionless function ^(e) is determined by a kinetic 
equation for the electron spin density matrix p nslS2 (z,p z ). 
It has been shown in Paper I, however, that a good ac- 
curacy is provided by a simpler kinetic equation for the 
density distribution function f ns = p nss in the "fixed spin" 



representation (Yakovlev 1984). 

In order to obtain the longitudinal transport coeffi- 
cients, it is sufficient to assume that E* and VT are 
collinear with B. The electron distribution function can 
be sought in the form 



fns = / + Zsign(p z )^r 
oe 



eE* 



e-n dT 
k B T ~dz 



Vns{e), (29) 



where / is an appropriate scale length, and 93(e) is a di- 
mensionless function to be determined from the kinetic 
equation. The latter is reduced to an algebraic system 
(Yakovlev |l984| ): 



a i 7 J->„' s '(e) [fns(e) - 7¥W(<0] = L 

ns — >-n' s' 



(30) 



Here, dyj^ n , s , is a dimensionless scattering rate of an elec- 
tron from a state with quantum numbers n and s into 
a state with quantum numbers n' and s', with changed 
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(7 = —1) or unchanged (7 = 1) direction of motion 
along B. The summation is performed over n' < fT. max (e), 
7 = ±1, and s' = ±1 for ri > 1 (but s' = -1 if n' = 0). 
Note that the terms with n' = n, s' = s, and 7 = +1 
naturally vanish. 

Since the scattering potential is written in the form 
(p0|), it is convenient to choose / = m c c 2 hL0 c /(2imiZ 2 e 4 ). 
Then (Yakovlev |l984 ) 



an 



ai 7 s L n , s ,(e) = [4(6+ 1) Wj-^^^Ce), 
where e = e/m c 2 , p„ = p n /m e c, 



7(7) 



[(?+ l) 2 +-fpnPn 



+4b 2 nn'Q 2 + 46 [(e + l) 2 + 7P«P«'] Vnri Q 3 , 

Snn',1 -1 = 2b [ n 'PnQl + n Pn< Qu] - 4& 7 ^/™« 7 PnPra' Q3, 

and S^J, _ a _ s , differ from S^,^/ by i nterch anging 



and Q 2 - Here the functions Qi (Yakovlev |l984|) are gener- 
alized to arbitrary scattering potential: 



Qi = I n -i, n i-i{u)4> (u)du, 
Jo 

/>oo 

Q 2 = / Z 2 Il ,(M)^ 2 (-u)dM, 
Jo 



(«)/„_!,„/_! (ti) 2 (u) du, 



where 

Inn' (u) 



1/2 



LI' 



(u) 



(33) 
(34) 
(35) 

(36) 



is a Laguerre function (Sokolov & Ternov 196S ; L™ (u) are 
the associated Laguerre polynomials - e.g., Abramowitz & 
Stegun 1972 ), (j){u) = 2\cj) q \/a 2 n , and (hq) 2 is set equal to 

(Pn ~ IPn') 2 + 2(fl/a m ) 2 U. 

3.4- Transport perpendicular to quantizing magnetic field 

Let us start with the case of large Hall parameter, tqlo & 
1. In this limit, as follows from Eq. (|l9|), 



l/wg and t xx w (u>Itq) l , 



(37) 



provided the magnetic field is non-quantizing. The case of 
quant izing field was con sidere d by Kaminker & Yakovlev 
( 1981 ) and Hernquist ( 1984 ). The expressions for the 
transport coefficients derived by these authors can be writ- 
ten in the form (fL8f) by defining, in analogy with Eq. (p7|), 



T yx = l/w g and t xx = (u 2 t±) \ 



(38) 



2 Equation (B2h reproduces Eq. (26) of Yakovlev ( |l984| ) cor- 
rected for a misprint. 



< 
&0 



(31) 5 



(32) 
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Fig. 2. Calculated (symbols) and fitted (solid lines) effec- 
tive longitudinal Coulomb logarithms A|| vs. dimension- 
less electron energy variable Eq. (^) for fully ionized iron 
(filled circles) and carbon (open circles) at B = 10 12 G and 
three values of temperature (lg T is given at the right end 
of each bunch of curves). Dot-dashed and dashed curves 
represent the field-free Coulomb logarithms Ao for Fe and 
C from Paper III. 



where 
N B (e)c 2 



Z 2 e i n i 



y( e ) 

e u 2 T ± (e) 27T?iV 2 1 >' 



*(e) 



y ^— 

*—f 2p n Pn' 
nn 7 

+ibVrm'Qj 



(f 2 + l + !PnPn-){Qi + Q2) 



(39) 



(40) 



and functions Qj- differ from Qi [Eqs. j3^)-(|35|)] by an 
additional factor u in each integrand. 

In weakly quantizing magnetic field, Tj_(e) oscillates 
around To(e); it can be replaced by 7o(e) in the non- 
quantizing limit. This allows us to interpolate between 
the regimes of large and moderate-to-low Hall parameters 
using the formulae: 



Tl 



1 + Kri) 2 



1 + (w g r ± ) 2 



(41) 



Equations ( j4l| ) relate the effective relaxation times t xx 
and r yx to the effective transverse electron collision fre- 



quency r^ 1 and correctly reproduce the known limits (19) 
(non-quantizing field, arbitrary Hall parameter) and (38) 
(arbitrary field, large Hall parameter). 

In addition, our interpolation ( |ll| ) of t xx eliminates 
the well known divergency, that arises from direct substi- 
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Fig. 3. Same as in Fig. ||but for B 
set of T. 



tution of Eq. ( J39| ) in the integrand of Eq. Jig ) because \P(e) 
turns to infinity at each Landau threshold as (y— n max )~ . 
Previously \&(e) was truncated at some level, estimated 
by a semi-qualitative analysis of physical processes that 
could, in principle, eliminate the divergency, were they 
included into the theory (Kaminker & Yakovlev 1981). 
One can show, however, that the derivation of expressions 
for the transverse transport coefficients [equivalent to our 
Eqs. (|3S ) — (p0|)] implied that <j s t± 3> 1. By correcting re- 
lations (38|) in case where w g rj_ is not very large, Eq. (^Tj) 
ensures finiteness of t xx , thus making a truncation unnec- 
essary. 



4. Fitting formulae for T||(e) and T±(e) 

Explicit expressions of Qi and Qj- for the function \<ff® | 2 
in the form (|2^) are given in the Appendix. Using them, 
we have performed extensive calculations of T\\ and tj_ 
from Eqs. (f|), @-@ and Q. The key parame- 
ters of the function \cf>° |, which enters these expres- 
sions, - the Debye-Waller parameter aow = w(2A;fo) ~ 

^2 



Fig. 4. Longitudinal (Ay , filled circles) and transverse 
(Aj_, open circles) effective Coulomb logarithms in iron 
plasma for B = 10 14 G and two values of T. The solid 
lines show the fit; the dot-dashed lines show Ao. 

Let us define longitudinal and transverse effective 
Coulomb logarithms Ay ,_i_(e) through the relations 



r±(e) 



Po v o 



M,(6) 

Af B (e) 47rniZ 2 e 4 A||(e)' 
P&o 



(42) 
(43) 



A/o(e) 4Trn i Z 2 e 4 A A _(e)' 

As seen from Eq. (pi]), the functions Ay and Aj_ turn into 
Ao if the magnetic field is non-quantizing. In the quantiz- 
ing fields, the ratios Ay/Ao and Aj_/Ao are fitted by the 
expressions 



Ay 

a7 



D 



2^2 



-L z x 



Po 
3x 2 



By/X + C- 



2n max + 1.5x 2 /(l + 2&) 2 



E 

0.07+ - 

5 



-1/2 



(44) 



Ax A A 

A« ^ 2 -2 



10 



u_2 (2fc.Fo/<7D) an( l the Coulomb screening parameter 
o-s = (?s/2A:fo) 2 , _ as well as the magnetic field param- 
eter b, varied independently from 10 -4 to 10 2 . For any 
value of b, the variable v [related to e via Eq. (^)] var- 
ied from to 25, taking on 5-10 values over each interval where x — p n /Vb = \J1{y 

Calculation is quite simple as long as v < 1, the func- a 
tions $(e) and \I/(e) being given by Eqs. ( A.8)-( A. 10). For 
v > 1, we have fitted the results of our numerical calcula- 
tions by analytic formulae. 



Po x" 

(B lnn max ) x^ 1 - (C + D lnn n 



)y/E , (45) 



c), 



30-15£;-(15-6£) vl 
30 - 10 E - (20 - 5 E) v 2 ' 



B = - — — 
2 2 



4 1 - 2w 2 /3' 



C = 



l-E + MSv 2 
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Fig. 5. Longitudinal (||) and transverse (_L) electrical (a) 
and thermal (k) conductivities in the outer neutron-star 
envelope composed of iron for B = 10 12 G and two val- 
ues of lgT (marked near the curves): comparison of the 
new results (solid lines) with the classical approximation 
[Eq. (|l9|), dashed lines] and with the old results for the 
longitudinal conductivities (dotted lines) from Paper II. 



Fig. 6. Same as in Fig. | but for B = 10 14 G and T 
10 8 K. 



D = 1 + 0.06- 



L' 2 



E 



1 - exp(-a DW ) 



A = 0.8(1 + L) + 0.2 L, B = (0.68 -1.3E) L 1/6 , 

C = 1.42 - E + £ 1/2 /3, D = fo.52 - Ej Z 1/4 , 

E = (10 + 5/6)~\ L = ln(l + a -1 ), L = aL, 

a = [V&T7 + (2 + 0.5 anwT 1 ] 2 ■ 

The leading terms at x — > are proportional to x^ 1 in 
Eq. ( p4| ) and x~ 2 in Eq. (Ejj), reproducing the asymptotic 
behaviour of the functions $ and "J. The accuracy of the 
fit was checked for v — n max > 0.01, which is quite suf- 
ficient for most applications. Equations (Q) and (45) fit 
our numerical results with a typical error of a few percent. 
A maximum error up to 40% occurs only at some extreme 
values of a s , adw, and v. In addition, we have compared 
the new fit ( ^i[ ) with the one in Paper II (in the particu- 
lar case of odw oo which corresponds to the screened 
Coulomb scattering potential considered in Paper II). On 
average, the new fit turned out to be more accurate than 
the old one. 

Figures |]-[| illustrate the accuracy of Eqs. (H) and 
( [45| ) for realistic parameters in the neutron-star envelopes. 
In Figs. [2] and [| the filled and empty circles represent 



A||(e) calculated from Eqs. (|28| ) and (^4) ; in Fig. [|we have 



additionally plotted Aj_(e), obtained using Eqs. ( p9|) and 
([43|). The parameters a s and odw have been calculated 
for fully ionized iron or carbon plasmas at various T from 
10 5 to 10 8 K and B from 10 12 to 10 15 G (indicated in the 
figures). At every point, the plasma density has been de- 
termined from the condition ep = e{y) using Eq. (0). In all 
figures, the fits ( [44] ) or (|4|) are drawn by solid lines. The 
dashed and dot-dashed lines represent the non-magnetic 
Coulomb logarithm Ao(e) given in Paper III. 

Since our fitting formulae depend analytically on the 
parameters of the effective screening function ([23]), they 
need not be changed in case future refinement of the the- 
ory will cause modification of these parameters. 

5. Numerical results for transport coefficients 

Figures ^| and ^ show electrical and thermal conductivities 
calculated with the effective relaxation times (e) given 
by Eqs. (|4l])-(^3|) at various p, T, and B appropriate for 
outer envelopes of th e neutr on stars. The use of the ana- 
lytic equations (|A^)- (|A.1QD for v < 1 and @, (|f) for 
v > 1 reduces numerical calculation to one-dimensional 
integration in Eq. (|l8|), which has been performed using a 
fast algorithm described in Sect. 5 of Paper II. 

The density range in every figures allows to see the 
strongly quantizing (below the first Landau threshold) and 
weakly quantizing regimes. The non-quantizing (classical) 
results are plotted by the dashed lines. 

Figure ^| shows the longitudinal and transverse con- 
ductivities in a neutron star envelope composed of iron 
for B = 10 12 G. The quantum oscillations around the 



A.Y. Potekhin: Electron conduction in magnetized neutron star envelopes 



795 




Fig. 7. Longitudinal conductivities for Z = 6 and T = 
10 6 K at B = 10 12 and 10 13 G. As in Figs. | and g, the 
solid, dotted, and dashed lines show the new, old, and 
classical results, respectively. 

classical values are more pronounced at lower tempera- 
tures. Figure (] illustrates the conductivities at stronger 
field, B = 10 14 G, which may be relevant to magnetars. 
The classical formulae correctly reproduce the large-scale 
trend of the curves and the reduction of transverse con- 
ductivities with respect to longitudinal ones. Nevertheless, 
deviations caused by the quantum oscillations are quite 
prominent, especially in the regime of strong quantization, 
where they may reach orders of magnitude. 

For comparison, we have plotted (by dotted lines) the 
longitudinal conductivities calculated using the formalism 
of Paper II. The temperature values in Fig. g have been 
deliberately chosen the same as in Fig. 5 of Paper II (10 6 
and 10 7 K). Note that for T = 10 6 K, the plasma is entirely 
within the solid crust. In this case, our old results agree 
nicely with the new ones. On the contrary, for T = 10 7 K 
in Fig. g as well as for T = 10 8 K in Fig. g, the displayed 
density range extends into both the solid crust and liq- 
uid ocean of the star. The new conductivities go smoothly 
across the phase transition, whereas the old ones exhibit 
large jumps and appear to be significantly overestimated 
just behind the ocean/crust interface. This is caused by 
an overestimated effect of the Debye- Waller factor in Pa- 
per II, now corrected by including multi-phonon processes. 

Figure g shows the longitudinal conductivities of car- 
bon plasma for B = 10 12 and 10 13 G, which may chance, 
e.g., in a neutron star with an accreted carbon shell. 
In this case, the bottom of the ocean lies slightly above 
p = 10 5 gem -3 . Once again, we observe significant discon- 



Fig. 8. Longitudinal (solid lines) and transverse (dashed 
lines) thermopower in units of fce/e, for two sets of indi- 
cated parameters. For comparison, the dotted lines show 
the results of Paper II. Vertical dot-dashed lines indicate 
the liquid/solid phase transition. 

tinuities of the "old" conductivities, which deviate from 
the new ones on both sides of the interface. The difference 
in the liquid phase is attributed to the modified ion struc- 
ture factor used for obtaining |0q ff | 2 in Paper III, instead 
of a simplified screened-Coulomb model in the previous 
work. 

In Fig. g we present the longitudinal and transverse 
components of the thermopower tensor for the sets of pa- 
rameters used in Figs, g, g The longitudinal thermopower 
from Paper II is shown by the dotted lines. Unlike the con- 
ductivities, the longitudinal thermopower did not possess 
considerable breaks at the phase transitions in the old the- 
ory. Nevertheless, one can observe that the "new" results 
differ from the "old" ones. As in the previous figures, this 
difference is noticeable in the vicinity of the freezing point. 

In the previous figures we have presented the longitu- 
dinal and transverse transport coefficients, which evince 
magnetic quantum oscillations around their classical val- 
ues. The off-diagonal (Hall) electrical and thermal con- 
ductivities do not exhibit such oscillations and practically 
coincide with their non-magnetic counterparts given by 
Eqs. (|l8|) and (|l9|). Unlike them, the Hall component of the 
thermopower, Q yx , does oscillate, as illustrated in Fig. g 
for iron plasma at B = 10 12 G and T = 10 6 , 10 7 , and 
10 s K. The oscillations are very sharp at T = 10 6 K, but 
they are completely smeared out at the highest temper- 
ature, T — 10 8 K, which is close to Tb in the present 
example. The difference of vertical scales in Figs, g and g 
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Fig. 9. Hall component of thermopower in units of k^/e, 
for iron plasma at B = 10 12 G and three values of tem- 
perature. 



reflects that Q yx is relatively small. Nevertheless, it may 
cause a variety of thermomagnetic effects in neutron star 
envelopes (Urpin et al. 



1986) 



6. Conclusions 

We have derived practical expressions for the electron 
transport coefficients in degenerate (T ^ Tp) layers of 
neutron-star envelopes with magnetic fields which may be 
quantizing. Generally these expressions require energy in- 
tegration, Eq. (|l8|). In the case of strongly degenerate elec- 
trons (T <C Tp) and not too close to the Landau thresholds 
[v— n max k-BT/hu; g ), even this numerical integration be- 
comes unnecessary, and the electrical conductivity reads 
<xy w (e 2 n e c 2 /cf) Uj (ep), where Tij is provided by our an- 
alytic formulae, while the thermal conductivity is given by 
the Wiedemann-Franz law, w (tt 2 k^T / 3e 2 ) (Jij . 

We have used an effective potential, which has been ob- 
tained in Paper III assuming that the conductivities are 
determined by the electron scattering off ions (off phonons 
in the cr ystalli ne phase). It is well known (e.g., Yakovlev 
& Urpin |198CD that this is the main mechanism regu- 
lating electron transport at T p ^ T ^ Tp. Other con- 
tributing mechanisms are the electron-electron scattering 
and scattering off lattice defects and impurities in the 
crystal. Corrections due to the impurity scattering can 
be introduced in a standard albeit approximate way by 
summation of relevant partial collisional frequencies, viz: 



1,-L 



unpi 



In the case of charged impuri- 
ties with charge number Z- lmp , occasionally embedded in a 
Coulomb lattice, the effective scattering potential is again 
given by Eqs. @ and (§|), by setting G = 1, e~ w ^ = 0, 
and qi = and replacing Z by | Z — Z- lmp | and n\ by n- lmp . 
Then Tm™l are given by our formulae with an obvious mod- 
ification of parameters. 

We expect that our new formulae for the conductivities 
will be useful, in particular, in calculations of neutron-star 
thermal structure and evolution. It would be especially 
interesting to apply these results to investigating thermal 



structure of magnetars. Up to now, a very simplified an- 
alytic model (Heyl & Hernquist 1998) has been used in 
this case, but the problem deserves a more thorough study 
since it may provide a clue to the origin of the anomalous 



X-ray pulsars (Heyl & Hernquist 1997). The formulae pre- 
sented here are almost as simple as those used by Heyl & 
Hernquist ( 1998| ), but they are accurate over a consider- 
ably broader range of plasma parameters. 

The computer code that implements the formulae de- 
rived in the present paper is freely available from the au- 
thor by electronic mail. 
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Appendix: calculation of auxiliary functions 



Let us adopt Eq. (|23j) and set G = 1 for simplicity. Then 
the dimensionless function 4> 2 {u) that enters Eqs. (^3|)- 
( [§5l) can be written as 

1 _ C -C(«+?) 

<t> {u) = . - V2 , (A.l) 

where ( = 2iu(a- 1 ) = 2(a m g D )- 2 u_ 2 (l + ,3/3), £ = 
\{a m /fi) 2 (p n - -ip n ') 2 , ito = £ + £ s , and £ s = \{a m q s ) 2 . 
Let us also define 



(o) f°° e -c(u+t) 

Qjnn'mi^ = J J «n' WW-jW ^ + 

Then, for j = or 1 , we have 



■du.(A.2) 
(A.3) 



Qt +J = Qfi'iK, 0) - u Q% n , 2 {u , 0) 



(«) 



a « 3 



OS»a(«o.O-«oQS°U(«o.O . (A-4) 



and Qi,Qj~ are obtained from Qi,Q 2 by replacing n — > 
n — 1 and n' — » n' — 1. 

For small Landau numbers (n, n' & 10), one can cal- 
culate Qjnn'm(£> using an explicit expression of the La- 



* jnn'm v 

guerre functions I nn '(u) in Eq. (A.2). Since I nn '(u) = 
(—1)" ~ n I n > n (u), we assume n' > n without any loss of 
generality. Then 



W«) = e^V"'""'/ 2 J2(-l) k Cn> nk i 



(A.5) 



k=0 



where c n i n k — y/n'\n\/[k\(n — k)\{n' — n + fc)!]. For j = 
or 1, 



2n~j min(n,/) 

°U(^c) =£(-!)< £ 

/— /c-max(0,/-n) 
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XC„>nkC n >,nd-k(n - 71 + l)\ Q^_„ + ;, , m (£, C), (A.6) 

S,0,i& = (! + 0""e^n+i(e + CO. 



,(0) 



0)2 (C, C) = (i + C) 1- "^^^ + CO - £n+i(£ + CO]; 

an exponential integral E n (x) — t~ ~ n e~ xt dt is easily 
calculated (Abramowitz & Stegun 1972). 

In case where n or n' is large, Eq. ( gjg ) is impracti- 
cal because of approximate cancellations of positive and 
negative terms. In this case, one can use the following rep- 
resentation^] 



} (0) 

■s jn'nm 



K.0 = E 



[n!(n')!(n~j)!(n'-j)!] 1/2 
fc! (n — j — fc)! (n' — j — fc)! 



" C ^e-^dx.(A.7) 



(m - 1)! (j + fc)! 7 C (1 + x ) n '+"-J+ 

Finally, let us consider an important particular case 
of n = n' = 0. For transport along and across magnetic 
field, the effective inter-collision times Tii(e) and Tj_(e) are 
related to the functions $(e) and \&(e) by Eqs. ( p8| ) and 
(p9|), respectively. From Eqs. and (E3), we obtain 



d>(e) = 



where 



*8 



2QH0 



vI/(e) = -[e 2 Q ± (0+g J 



(A.8) 



Q"(0 = u - 1 (l-e-«)-e Uo Bi(uo) 



+(1 + C) c" 0+ces Si(u + C«o), 



(A.9) 



Q (0 = (l+u )e t,0 Ei(uo)-l + ' 



-(1 + u + Cuo)e 



£i(u + C"o), (A. 10) 



£ = 2pq/&, and the function _Ei is readily given by poly- 
nomial approximations (Abramowitz & Stegun 1972). 
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